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Reaction  Kinetics  /  Unimolecular  Reactions  /Ion  Fragmentations 


The  distribution  of  translational,  rotational,  and  vibrational  energies  in  the  fragments 
(benzylium  ions  and  ethyl  radicals)  of  the  dissociation  of  »-propylbenzene  cations  has 
been  determined  by  statistical  adiabatic  channel  model/classical  trajectory  (SACM/CT) 
calculations.  The  reaction  was  treated  by  CT  calculations  of  capture  processes  for 
transitional  modes,  starting  with  specified  fragment  energies.  A  short-range  valence/long- 
range  ion-induced  dipole  potential  model  for  the  transitional  modes  was  employed. 
The  derived  distributions  approach  the  results  from  phase  space  theory  (PST)  at  small 
energies  and  angular  momenta.  At  larger  energies  and  angular  momenta,  the  shapes 
of  the  distribution  functions  remain  similar  to  those  from  PST;  however,  the  average 
translational,  rotational,  and  vibrational  energies  of  the  fragments  increasingly  differ  from 
PST  predictions.  The  present  results  are  consistent  with  separate  SACM/CT  calculations 
on  the  same  potential  of  specific  rate  constants  k(E,  J)  and  thermally  averaged  rate 
constants  kx(T)  of  the  dissociation/recombination  reaction. 


1.  Introduction 

Translational  product  energy  distributions  (TPEDs)  constitute  important  ob¬ 
servables  characterizing  the  dynamics  of  unimolecular  dissociation  processes. 

*  Corresponding  author.  E-mail:  shoff@gwdg.de 


i 

2006011/  465 


700 


J.  Troe  et  al. 


Some  information  exists  for  dissociations  of  neutral  molecules,  but  TPEDs 
have  most  frequently  been  measured  for  the  fragmentation  of  molecular 
ions.  The  available  experimental  techniques  have  well  been  described  in 
the  monograph  by  T.  Baer  and  W.  L.  Hase  [1].  Particularly  important  tech¬ 
niques  that  measure  TPEDs  are  PEPICO  (photoelectron  photoion  coincidence) 
and  MIKES  (mass  analyzed  ion  kinetic  energy  spectroscopy).  The  measured 
TPEDs  have  been  analysed  by  a  variety  of  models,  which  also  have  been 
summarized  in  [1],  The  most  useful  schemes  have  been  orbiting  transition 
state/phase  space  theory  (OTS/PST)  such  as  elaborated  by  Chesnavich  and 
Bowers  (see  e.g.  [2-5]),  phase  space  theory  involving  the  concept  of  product 
temperatures  (PST,  T*)  by  Klots  (see  e.g.  [6-8]),  variational  flexible  transition 
state  RRKM  theory  with  exit  channel  coupling  (VRRKM/ECC)  by  Wardlaw, 
Marcus,  Klippenstein,  et  al.  (see  e.g.  [9-11]),  the  separate  statistical  ensem¬ 
bles  (SSE)  model  by  Wittig  et  al.  [12],  and  the  statistical  adiabatic  channel 
model  (SACM)  [13].  Comparison  with  experimental  TPEDs  often  showed 
good  agreement  with  OTS/PST  calculations.  Dissociations  with  large  barriers 
in  the  reverse  direction  gave  rise  to  enhanced  kinetic  energy  release  (KER). 

The  success  of  OTS/PST  may  appear  troubling  for  the  following  rea¬ 
sons.  While  TPEDs  appear  to  be  modeled  adequately  by  loose  activated  com¬ 
plex  PST,  specific  rate  constants  k{E)  of  ion  fragmentation  traditionally  are 
represented  by  rigid  activated  complex  RRKM  theory  with  oscillator-type 
transition  states.  k(E)  as  well  as  the  corresponding  thermally  averaged  high 
pressure  dissociation  rate  constants  kx{T)  are  often  far  below  OTS/PST  cal¬ 
culations.  Thus,  there  seems  to  be  an  inconsistency  which  deserves  further 
investigation.  The  present  article  addresses  this  subject  for  the  fragmentation  of 
n-propylbenzene  cations  as  a  representative  example  of  a  simple  bond  fission 
without  reverse  barrier.  The  reaction  predominantly  proceeds  via  a  C-C  bond 
fragmentation, 

h-C9H12+-*C7H7+  +  C2H5.  (1.1) 

Specific  rate  constants  k(E )  for  this  reaction  have  been  measured  by  Kim 
et  al.  using  laser  photodissociation  and  MIKES  [14]  and  we  have  recently 
measured  thermal  dissociation  rate  constants  [15],  k(E)  was  analysed  in  [14] 
by  rigid  activated  complex  RRKM  theory,  whereas  TPEDs  were  represented 
by  OTS/PST.  In  contrast  to  this,  we  modeled  k(E,  J )  by  SACM/CT  (sta¬ 
tistical  adiabatic  channel  model/classical  trajectories)  calculations  [15, 16]. 
While  it  appears  difficult  to  combine  RRKM  and  OTS/PST  into  one  inter¬ 
nally  consistent  scheme,  this  is  not  the  case  with  SACM/CT.  This  approach 
allows  one  to  generate  k(E,  J )  and  TPEDs  as  a  function  of  energy  E  and  an¬ 
gular  momentum  (quantum  number  J)  in  a  consistent  manner.  The  present 
work  describing  SACM/CT  calculations  of  TPEDs  for  reaction  (1.1)  comple¬ 
ments  [15]  and  [16]  where  k(E,  J )  and  kx(T)  were  modeled  using  the  same 
potential. 
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2.  Calculational  method 

2.1  SACM/CT  concept 

Since  the  concept  of  the  statistical  adiabatic  channel  model  often  is  misrepre¬ 
sented  (even  in  the  monograph  [1]),  the  principles  of  our  approach  are  briefly 
summarized  in  the  following:  (i)  The  internal  modes  of  the  reactant  and  the 
fragments  are  compared  in  order  to  identify  conserved  and  transitional  (dis¬ 
appearing)  modes.  The  conserved  modes,  to  a  first  approximation,  are  identi¬ 
fied  with  the  internal  modes  of  the  fragments.  In  dissociation,  the  conserved 
modes  are  assumed  to  stay  adiabatically  in  their  respective  quantum  states 
during  the  final  separation  of  the  fragments  (or  during  the  initial  approach  of 
the  fragments  in  the  reverse  association).  In  other  words,  our  approach  as¬ 
sumes  vibrational  adiabaticity  in  the  conserved  modes,  (ii)  The  dynamics  of 
the  dissociation  is  divided  into  two  stages.  During  the  first  stage,  intramolec¬ 
ular  vibrational  redistribution  (IVR)  or  other  processes  (such  as  collisions) 
are  assumed  to  establish  a  statistical  distribution  of  reactant  states.  During  the 
second  stage,  the  fragments  must  pass  dynamical  bottle-necks  created  by  the 
anisotropy  of  the  potential  and/or  by  centrifugal  effects  before  separation  oc¬ 
curs.  Not  all  trajectories  are  successful  at  passing  the  bottle-neck.  The  reverse 
association  starts  from  a  statistical  distribution  of  reactants  and  the  combin¬ 
ing  fragments  must  pass  the  same  dynamical  bottle-necks  for  association  to 
occur.  In  the  second  stage  then  capture  is  terminated  by  IVR  or  by  collisions. 
The  assumptions  of  vibrational  adiabaticity  and  of  statistical  distributions  at  the 
beginning  of  the  final  separation  of  the  dissociation  fragments  provide  a  con¬ 
siderable  simplification  of  the  treatment  of  the  complete  dynamics.  However, 
these  assumptions  are  not  necessarily  valid  and  need  to  be  justified,  (iii)  The 
dynamics  of  the  transitional  modes  during  the  second  stage  of  dissociation  (or 
the  reverse  first  stage  of  association)  is  treated  explicitly  on  a  potential  energy 
surface  (PES)  of  reduced  dimensionality,  i.e.  on  a  PES  with  frozen  conserved 
modes.  The  dynamics  of  transitional  modes  may  be  rotationally  adiabatic  or 
nonadiabatic,  in  the  limit  being  “sudden”,  see  below.  The  extent  of  the  nonadia- 
baticity  is  governed  by  the  effective  mass  M  of  the  system  or  the  corresponding 
Massey  parameter  £  =  a/2 M  which  corresponds  to  the  ratio  of  the  effective  ro¬ 
tational  period  of  the  transitional  modes  and  the  collision  time  on  the  reduced 
dimensionality  PES.  In  the  adiabatic  limit  (i.e.,  in  the  limit  £  »  1),  the  motion 
proceeds  on  one-dimensional  adiabatic  channel  potential  curves  which  can  be 
constructed  by  determining  the  eigenvalues  of  the  transitional  modes  at  frozen 
interfragment  center-of-mass  distances.  In  the  earliest  version  of  SACM  [17], 
these  eigenvalues  were  approximated  by  exponential  interpolations  between 
reactant  and  fragment  eigenvalues  assuming  a  non-crossing  rule.  In  later  ver¬ 
sions  of  SACM  (see  e.g.  [18-23]),  eigenvalues  were  calculated  accurately  for 
simple  model  potentials.  There  are  avoided  crossings  of  the  adiabatic  channel 
potentials,  few  for  atom  +  linear  capture  and  numerous  for  linear  +  linear  cap¬ 
ture.  Nonadiabatic  transitions  occur  at  such  crossings  or  over  broader  ranges  of 
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distances,  Coriolis  coupling  being  of  particular  importance  [24],  Nonadiabatic 
dynamics  (in  the  limit  being  “sudden”,  i.e.,  £  1)  most  conveniently  is  treated 

by  classical  trajectory  (CT)  calculations.  It  was  shown  [25]  that  adiabatic  chan¬ 
nel  and  CT  treatments  in  the  range  of  classical  and  adiabatic  dynamics  lead 
to  identical  results.  Therefore,  SACM/CT  appears  to  be  the  method  of  choice. 
Refined  versions  of  variational  transition  state  theory  that  include  dynamical 
corrections  [26]  and  the  SACM/CT  outlined  here  appear  to  yield  similar  results 
if  they  are  applied  to  the  same  PES.  However,  in  both  cases  CT  calculations  are 
essential  for  an  accurate  characterization  of  the  transitional  mode  dynamics. 

Specific  rate  constants  k(E,  J )  within  the  SACM/CT  concept  are  repre¬ 
sented  by  the  well-known  relationship  from  statistical  unimolecular  theory  [1] 

k(E,  J )  =  W(E,  J)/hp{E,  J)  (2.1) 

where  p(E,  J)  denotes  the  rovibrational  density  of  states  of  the  reactant 
molecule.  The  contribution  of  the  transitional  modes  to  W(E,  J)  in  SACM/CT 
is  given  by  the  number  of  quantum  states  W0(E,  J )  of  the  separated  fragments 
at  a  given  E  and  J  multiplied  with  the  capture  probability  w(E,  J)  such  as 
derived  by  CT  calculations.  In  the  limit  of  adiabatic  dynamics,  this  product 
corresponds  to  the  number  of  adiabatic  channels  which  are  open  in  associa¬ 
tion  or  in  the  reverse  dissociation.  In  the  case  of  nonadiabatic  dynamics,  axially 
nonadiabatic  versions  of  SACM  [24]  or  CT  [25]  calculations  are  used  for  the 
determination.  Compared  to  a  variational  treatment  applied  to  phase  space  vol¬ 
umes  such  as  used  in  VTST,  the  SACM/CT  approach  has  the  advantage  of 
automatically  including  dynamical  corrections.  The  total  W(E,  J)  finally  is  ob¬ 
tained  by  convolution  of  the  contribution  from  the  transitional  modes  with  the 
contribution  from  the  quantized  conserved  modes  such  as  this  is  usual  in  statis¬ 
tical  unimolecular  rate  theory  [1], 

For  the  present  purpose,  Eq.  (2.1)  is  further  specified  by  selecting  trajecto¬ 
ries  in  CT  capture  calculations  with  a  given  property  of  interest,  e.g.  a  specific 
product  translational  or  rotational  energy.  Collecting  those  initial  states  from 
the  phase  space  of  the  separated  fragments,  e.g.,  with  given  translational  ener¬ 
gies  E„  which  lead  to  capture,  gives  W(E,  J,  Ea).  The  derived  TPEDs  then  are 
given  either  in  cumulative,  integral,  form  as 


Eu  E 

P  (E,  J,  Ea)  —  J  W  (E,  J,  E'tt)  d E'a /  J  W(E,  J,  E’tt)  d E’tt  (2.2) 

o  o 

or  in  differential  form  by  dP(E,  J,  Ett)/dEa.  For  calculational  purposes  and 
representational  reasons  in  the  present  work  we  prefer  the  integral  form  of 
Eq.  (2.2)  over  the  differential  form.  Vibrational  and  rotational  energy  distri¬ 
butions  can  be  determined  in  an  analogous  manner.  For  details  of  our  CT 
calculations,  see  [25], 
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The  described  SACM/CT  calculation  of  TPEDs  is  consistent  with  the 
SACM/CT  determination  of  k(E,  J)  and  of  kx{ T).  Our  concept  of  determin¬ 
ing  TPEDs  from  CT  capture  calculations  is  not  completely  new.  It  corresponds 
to  the  rudimentary  treatment  for  a  simple  triatomic  system  by  Hamilton  and 
Brumer  [27].  It  should  finally  be  emphasized  that  in  applying  CT  calculations 
to  the  transitional  mode  dynamics,  zeropoint  energy  problems  are  negligible 
for  most  conditions  of  interest.  However,  at  ultralow  temperatures  or  energies 
very  close  to  threshold,  the  corresponding  quantum  effects  could  also  be  ac¬ 
counted  for  in  the  usual  way  by  SACM  or  by  axially  nonadiabatic  SACM  such 
as  discussed  in  [24]. 


2.2  Potential  energy  surface 

The  described  SACM/CT  concept  can  be  applied  to  ab  initio  PESs,  whenever 
these  are  available  [28],  or  to  modeled  PESs.  For  the  CT  part  of  the  treatment, 
reduced  dimensionality  PESs  of  the  transitional  modes  are  required.  Ab  initio 
PESs  generally  are  not  yet  available  or  accurate  enough  in  die  range  of  in¬ 
terfragment  center-of-mass  distances  r  where  the  dynamical  bottle-necks  are 
located  (see  e.g.  the  N4+  potential  considered  in  [29]  and  [30]).  Therefore, 
one  often  has  to  live  with  modeled,  hopefully  realistic,  PESs.  In  the  present 
case,  we  had  to  construct  a  model  potential  which  reproduces  the  essential  fea¬ 
tures  of  C9H12+  and  of  the  fragments  C7H7+  4-  C2H5  and  which  characterizes 
their  interaction  in  a  simplified  manner.  This  required  at  least  one  fit  parame¬ 
ter  which  was  adjusted  by  comparison  with  experimental  k(E).  After  that,  no 
further  adjustments  were  made.  Our  PES  was  described  in  detail  in  [16]  and  is 
only  briefly  summarized  in  the  following. 

We  separate  the  PES  of  the  transitional  modes  into  the  sum  of  a  radial  and 
an  angular  part,  V(r)  and  V(r,  angles)  respectively.  The  radial  part  is  character¬ 
ized  by  a  switching  between  a  short-range  valence  potential  VSR(r)  of  Morse 
type  and  a  long-range  modified  electrostatic  (in  the  present  case  ion-induced 
dipole)  potential  FLM(r)  where 

VSR(r)  =  E0{  exp[— 2/5 (r  -  re)]  -  2  exp[-/3(r  -  r.)]}  (2.3) 


and 


VLR(r) 


aq2/2 

(r-  ArB)4+aq2/2D 


(2.4) 


with  A rB  =  re  —  rBe.  rBe  represents  the  equilibrium  value  of  the  length  of  the 
forming  or  breaking  bond  and  rc  is  the  equilibrium  value  of  the  c.o.m.-distance 
between  the  fragments.  At  large  r,  i.e.  r  »  ArB  and  r4  »  aq2/2D,  VLR{r)  ap¬ 
proaches  the  ion-induced  dipole  potential  —aq2/2r4.  The  switching  between 
VSR(r)  and  ELR(r)  is  performed  by  the  expression 


V(r)  =  FSR(r)[l  -kF}+  VLR(r)F 


(2.5) 
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with 

F=  l/{l+exp[-2£(r-rsw)]}  (2.6) 

and 


aq 2 

E0(rsv/  -  A rB)4  ' 


(2.7) 


rsw  denotes  that  value  of  r  where  VSR(r)  and  VLR(r)  from  Eqs.  (2.3)  and 
(2.4)  intersect.  The  relevant  parameters  for  reaction  (1.1)  were  chosen  as 
Eo/hc  =  13950cm-',  p  =  5.13  A"1,  re  =  4.74  A,  o^H;)  =  4.27  A3,  rBe  = 
2.42  A,  rsw  =  5. 14  A,  see  [16], 

The  angular  potential  is  represented  by  an  anisotropy  of  dipole-dipole  type, 
with  an  amplitude  which  exponentially  decreases  with  increasing  r,  through 

V(r,  angles)  =  C  exp  [  -  /I  (r  —  r  e)  ]  { [2  +  </,  </2  -  3  (rf,  n)  (</2  n)  ]  /4 } .  (2.8) 

d\  and  d2  are  unit  vectors  along  the  axes  of  the  fragments,  n  is  the  unit 
vector  along  the  line  connecting  the  c.o.m.s.  The  amplitude  of  the  anisotropy 
in  principle  corresponds  to  that  of  the  valence  potential.  In  practice,  however, 
we  treat  this  as  a  fit  parameter  to  reproduce  measured  values  of  k(E).  In  the 
present  case,  for  illustrative  purposes,  C/E0  ~  20  was  chosen  like  in  [16],  The 
value  of  C  is  of  particular  importance  since  it  defines  the  deviation  from  PST 
which  is  characterized  by  C  —  0;  the  particular  shape  of  the  angular  depen¬ 
dence  in  Eq.  (2.8)  is  only  of  much  smaller  relevance.  E.g.,  in  [16]  we  also 
studied  anisotropies  of  ion-permanent  dipole  type  which  led  to  similar  results, 
however,  with  a  different  anisotropy  amplitude  factor. 

In  order  to  simplify  our  CT  calculations  in  the  extensive  calculations 
described  below,  both  C7H7+  and  C2H5  were  treated  as  linear  rotors  with 
the  effective  rotational  constants  BJhc  =  (0.18  x  0.094  x  0.062) 1/3  cm-1  = 
0.102  cm-'  and  B2/hc  =  (3.51  x  0.753  x  0.699) 1/3  cm-1  =  1.23  cm-1  respec¬ 
tively.  In  order  to  classify  the  capture  dynamics  of  the  transitional  modes  with 
respect  to  the  degree  of  rotational  adiabaticity,  we  consider  an  average  Massey 
parameter  which  is  given  by  [16] 

f  ^  (2 iiB)l<2/hp  (2.9) 

with  S  «  (fi,  +  B2)/2.  This  Massey  parameter  is  equal  to  f  ~  0.18.  Thus,  the 
dynamics  of  the  transitional  rotational  modes  is  not  adiabatic  but  on  the  way  to 
sudden  dynamics  (f  <$C  1).  On  the  other  hand,  the  dynamics  of  the  conserved  vi¬ 
brational  modes  is  adiabatic  (£  1).  This  is  found  by  replacing  B  in  Eq.  (2.9) 

by  vibrational  quanta. 
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3.  CT  calculations  of  product  energy  distributions 

3.1  Centrifugal  barriers 

Translational  product  energy  distributions  are  markedly  influenced  by  angu¬ 
lar  momentum  conservation  constraints  as  well  as  by  the  centrifugal  barriers 
Eq(  J)  arising  from  the  PES.  In  the  simplest  version  of  phase  space  theory, 
E0(J)  would  simply  be  put  equal  to  zero.  Larger  kinetic  energy  release  is 
obtained  in  OTS/PST  where  E0(J )  is  derived  from  the  electrostatic  long- 
range  potential.  In  SACM/CT  we  go  one  step  further  by  determining  E0(J) 
from  the  complete  PES  except  that  the  anisotropy  is  artificially  switched  off. 
A  radial  potential  is  employed  which  accounts  for  the  complete  short-range 
valence/long-range  electrostatic  potential.  This  contrasts  to  the  artificial  lim¬ 
itation  of  only  a  long-range  potential  used  in  OTS/PST.  In  other  words,  we 
calculate  E0(  J)  from  the  complete  minimum  energy  path  of  the  PES. 

E0(J)  was  either  obtained  by  CT  calculations  or  by  determining  the 
maxima  of  the  effective  potentials  V(r)  +  L(L  +  l)h2/2/ir2  with  V(r)  from 
Eq.  (2.5)  ( L  —  orbital  angular  momentum  quantum  number);  the  approxima¬ 
tion  J  «  L  in  the  present  case  can  be  used  for  most  conditions,  see  [16].  The 
results  from  both  methods  can  be  represented  by 

E0(J)/hc~5.90x  10~7  cm-1  {[./(./+ l)]2/[l +9.73  x  10_3J107S 

+  5.21  x  10_6J2 166]} .  (3.1) 

Eq.  (3.1)  illustrates  the  differences  between  our  PST  and  OTS/PST  for  which 
E0(J)  would  be  given  by 

EoW^lJiJ+Dhi/lvf/laq2,  (3.2) 

i.e.  from  Eq.  (3.1)  with  the  denominator  put  equal  to  unity.  One  notices  that 
the  short-range  valence  part  of  the  potential  at  large  J  reduces  E0(J)  in  com¬ 
parison  to  that  found  from  only  the  long-range  potential.  As  a  consequence  at 
high  temperatures,  the  present  k™J(T)  slightly  exceeds  the  Langevin  capture 
rate  constant,  kL  =  2nq(a/i.i)'/2,  which  would  be  the  result  from  OTS/PST, 
see  [16]. 


3.2  Specific  capture  probalitities  w(E,  J)  and  rigidity  factors  frig(E,  J) 

Product  energy  distributions  are  most  conveniently  interpreted  by  looking  first 
at  overall  capture  probabilities,  w(E,  J),  for  which  the  product  energies  are 
not  yet  specified.  For  this  reason,  in  Fig.  la-c  we  present  capture  probabili¬ 
ties  w(E,  J )  and  PST  capture  probabilities  wVST(E,  J )  for  J  =  0, 100,  and  200, 
respectively.  In  Fig.  2a-c  specific  rigidity  factors, 

fnS  =  w(E,  J)/wTST(E,  J), 


(3.3) 
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Fig.  1.  Capture  probabilities  w(E,J)  (full  lines)  and  capture  probabilities  wPST(E,  J) 
from  phase  space  theory  (dashed  lines)  for  the  dissociation  of  n-propylbenzene  cations. 
SACM/CT  calculations  from  this  work  for  J  =  0  (Fig.  la,  E0(J)  =  0),  J  —  100  (Fig.  lb, 
£o(-f)/^c  =  23.8cm_1),  and  .7  =  200  (Fig.  lc,  Efs{J)/hc  =  2\Acmr'). 
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are  plotted  as  a  function  of  the  energy  E  and  at  the  same  fixed  angular  momenta 
(J  =  0, 100,  and  200).  PST  where  centrifugal  barriers  are  neglected  would  give 
wTSJ(E,  J)  =  1,  such  that  W{E,  J)  =  W0(E,  J)  would  correspond  to  the  well- 
known  statistical  expressions  given  in  [5, 31-33]  (before  convolution  with  the 
conserved  modes).  wPST(E,  J)  in  the  present  case  can  well  be  approximated  by 

wrST(E, ./)  ~  (1  —  E0(J)/E)2  (3.4) 

which  corresponds  to  the  large  J-approximation  ( J  «  L).  According  to  the  dif¬ 
ferences  in  E0{J)  such  as  quantified  by  Eq.  (3.1),  OTS/PST  and  our  present 
PST  through  Eq.  (3.4)  have  slightly  different  values  of  wPST(E,  J).  Fig.  la-c 
illustrate  the  effects  of  E0(J)  on  wfST(E,  J).  The  decrease  of  w(E,  J)  from 
ujpst(£,  J)  arises  from  the  anisotropy  of  the  potential  such  as  illustrated  in 
Fig.  la-c.  The  observed  effects  are  easily  interpreted.  At  low  energies  not 
far  above  E0(J),  the  dynamical  bottle-necks  are  located  at  large  c.o.m.  dis¬ 
tances  where  the  PES  is  not  very  anisotropic.  With  increasing  energies  and 
increasing  J,  these  bottle-necks  move  to  smaller  values  of  r  where  the  in¬ 
creasing  anisotropy  enhances  the  constraints  and  w(E,  J )  increasingly  falls 
below  wVST(E,  J).  This  effect  is  seen  particularly  clearly  in  Fig.  2a-c.  One  may 
expect  that  the  corresponding  “rigidity  effects”  also  become  apparent  in  the 
TPEDs  described  in  the  following. 

3.3  Translational  product  energy  distributions 

Fig.  3a-c  compare  cumulative  TPEDs  defined  by  Eq.  (2.2)  for  the  anisotropic 
model  PES  of  reaction  (1.1)  with  the  corresponding  TPEDs  for  the  artificially 
isotropic  model  PES  characterizing  our  PST.  Our  CT  results,  for  fixed  J,  are 
plotted  as  a  function  of  the  product  translational  energy  divided  by  the  total  en¬ 
ergy  £total.  It  should  be  noted  that,  in  our  work,  the  zeropoint  of  the  energy  scale 
is  put  at  the  lowest  possible  rovibrational  energy  of  the  separated  fragments.  In 
order  to  emphasize  the  translational,  rotational,  and  vibrational  contributions 
to  the  energy,  the  energy  E  is  also  called  Etoai.  Near  to  E  =  0  where  vibra¬ 
tions  cannot  be  excited,  E  =  Etotal  —  E„  +  Emt  such  that  P(E,  J,  Emt)  =  1  - 
P(E,  J,  Ea).  For  larger  Evib,  the  relation  is  more  complicated  but  can  be  recon¬ 
structed  from  Figs.  3-5.  Vibrations  are  always  treated  quantized  while  rotations 
and  translations  are  classical. 

A  number  of  observations  are  worth  mentioning:  (i)  For  small  values  of 
Etotai,  the  TPEDs  do  not  differ  too  much  from  the  distributions  from  PST. 
However,  with  increasing  Etoal,  the  TPEDs  are  increasingly  moving  to  the 
left,  i.e.  there  is  increasingly  more  KER  than  in  PST.  (ii)  For  J  =  0,  the 
cumulative  TPEDs  can  be  representated  well  by  1  -  exp(— E'tr/(£'tr)).  This  cor¬ 
responds  to  exponential  differential  TPEDs,  exp (-Ea/{Ea}).  With  increasing 
J  and  small  £tota ,,  however,  this  changes  into  differential  TPEDs  of  the  type 
^"expC-^tr/ (£)r)).  (iii)  With  increasing  2stotal,  the  average  KER  (Ea)  rises  in¬ 
creasingly  above  ( E „)  from  PST.  Our  PST  results  correspond  to  the  typical 
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Fig.  3.  Cumulative  translational  product  energy  distributions  (Eq.  (2.2))  in  the  dissoci¬ 
ation  of  n-propylbenzene  cations.  SACM/CT  calculations  from  this  work  for  the  full 
anisotropic  potential  (full  lines)  and  for  the  artificially  isotropic  potential  (PST,  dashed 
lined).  Curves  from  left  to  right:  Em\/hc  cm-1  =  10000,  3160,  1000,  316,  100,  31.6,  10, 
3.16,  1,  and  0.316  for  ^  =  0  (Fig.  3a),  E^/hc  cm'1  =  10024,  5480,  2996,  1638,  896, 
490,  268,  146,  80,  43.7  for  ./  =  100  (Fig.  3b),  and  EMA/hc  cm"1  =  10214,  6939,  3203, 
2176,  1479,  1005,  682,  464,  315  for  J  =  200  (Fig.  3c). 
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behaviour  of  PST  such  as  described,  e.g.  in  [1-8].  The  deviations  from  the 
TPEDs  of  PST  are  caused  by  the  anisotropy  of  the  PES. 

We  do  not  intend  at  this  stage  to  provide  further  quantitative  analysis 
of  the  TPEDs  of  Fig.  3a-c.  However,  we  note  that  increasing  deviations  of 
TPEDs  from  OTS/PST  with  increasing  energy  have  also  been  observed  exper¬ 
imentally,  e.g.  in  the  dissociation  CH3I+  — »■  CH3+  +  I  [35,36].  Exponential 
differential  TPEDs  have  also  been  observed  experimentally  [35-38],  However, 
our  calculations  show  that  one  must  expect  deviations  from  this  behaviour  with 
increasing  J  and,  in  particular,  if  the  fragments  both  are  not  monatomic.  It  has 
been  suggested  that  ( Ett )  should  be  proportional  to  £total  [39]  (see  the  represen¬ 
tation  of  results  for  CH3COCH3+  CH3CO+  +  CH3  in  [40]).  Fig.  3a  shows 
that  there  is  a  range  at  small  Eloa]  (and  small  J )  where  such  a  behaviour  is  also 
observed  in  our  calculations.  However,  with  increasing  Eloa]  both  in  PST  and 
in  the  general  PES  the  ratio  (Ea) /Et0 tai  obviously  decreases,  finally  reaching 
a  stationary  value  of  ( Ea )  independent  of  Etaal. 


3.4  Rotational  and  vibrational  product  energy  distributions 

Our  CT  calculations  also  allowed  for  the  specification  of  product  rotational 
and  vibrational  energy  distributions.  We  illustrate  the  corresponding  cumula¬ 
tive  distributions  in  Figs.  4  and  5,  again  as  a  function  of  the  energies  relative 
to  Emat  and  as  a  function  of  angular  momentum.  The  quantum  nature  of  the 
distributions  after  convolution  with  the  conserved  modes  becomes  quite  appar¬ 
ent,  particularly  in  the  vibrational  energy  distributions,  see  Fig.  5a-c,  but  also 
in  the  rotational  distributions  of  Fig.  4a-c  where  the  most  pronounced  “irreg¬ 
ularities”  arise  at  an  energy  where  the  first  vibrational  levels  of  the  fragments 
can  be  excited. 

We  do  not  intend  to  further  analyze  the  fine  details  of  the  distributions  here 
but  we  emphasize  that  there  are  ranges  where  the  distributions  markedly  dif¬ 
fer  from  the  PST  distributions.  The  differences  are  not  so  much  in  the  shape 
of  the  distribution  functions  but  more  in  the  average  energies  {Em)/Etoat  and 
(£vib)/£totai-  The  specific  properties  of  the  vibrational  and  rotational  energy 
distributions  depend  on  the  specific  molecular  parameters  of  the  considered  re¬ 
action  systems.  Not  all  details  of  Figs.  3-5,  therefore,  may  be  expected  to  be 
found  in  other  reaction  systems. 


4.  Conclusions 

The  present  SACM/CT  calculations  of  product  energy  distributions  in  the  frag¬ 
mentation  of  «-propylbenzene  cations  have  shown  that  PST  provides  more  or 
less  realistic  shapes  of  the  distribution  functions.  However,  the  average  trans¬ 
lational,  rotational,  or  vibrational  energies  of  the  fragments  are  close  to  the 
values  from  PST  only  for  small  J  and  small  energies  E.  At  larger  E  and  J,  one 
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Fig.  5.  As  Fig.  3  for  vibrational  product  energy  distributions,  however,  with  curves  from 
right  to  left  for  the  indicated  values  of  Eaul . 
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has  to  expect  increasing  differences  from  PST.  This  behaviour  is  due  to  the  shift 
of  the  effective  bottle-necks  of  the  reaction  from  the  isotropic  long  range  (ion- 
induced  dipole)  part  of  the  potential  into  the  increasingly  anisotropic  region  of 
the  short-range  valence  potential.  The  corresponding  behaviour  has  also  been 
observed  in  our  previous  SACM/CT  modellings  of  the  specific  rate  constants 
k(E,  J)  and  the  thermally  averaged  dissociation  rate  constants  kx(T)  of  the 
reaction.  More  studies  of  product  energy  distributions  by  the  present  type  of 
SACM/CT  capture  calculations  for  other  dissociation  processes  appear  desir¬ 
able. 
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